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Weak interactions between quarks give rise to hadronic parity violation which can be observed in 
nuclear and few-nucleon systems. We study the QCD renormalization of the isotensor component of 
parity violation at next-to-leading order accuracy. The renormalization group is employed to evolve 
the interaction down to hadronic scales. As the results are renormalization scheme dependent, 
we compare various schemes, including 't Hooft-Veltman dimensional regularization, and several 
regularization independent-momentum subtraction schemes. 
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I. INTRODUCTION 

Measurement of parity violation in nuclear and few- 
nucleon reactions allows one to probe the flavor- 
conserving hadronic weak interaction. Experimentally 
hadronic parity violation has been observed in a variety 
of reactions, see [l|, 0] for reviews. The interpretation 
of experimental data has been carried out using various 
models of the parity-violating nuclear force, the most 
popular of which is the meson-exchange model of Des- 
planques, Donoghue, and Holstein In this model, the 
parity- violating nuclear force is parametrized in terms of 
several parity-violating couplings between nucleons and 
mesons, which then give rise to parity violation in nu- 
clei. Analyzed within this framework, constraints on 
parity-violating couplings coming from different exper- 
iments are not consistent. Such discrepancies could arise 
from several sources: uncertainty entering nuclear struc- 
ture computations, model assumptions about the parity- 
violating nuclear force, dynamical effects due to the non- 
perturbative nature of QCD, etc. Resolving this situa- 
tion and connecting parity violation observed in nuclear 
reactions to parameters in the Standard Model is an am- 
bitious undertaking. 

Great progress has been made recently in describ- 
ing parity-violating few-nucleon reactions in a model- 
independent fashion [34l0j. Ever advancing experiments 
in few-body systems, moreover, seek to better constrain 
the parity- violating nuclear force. Bounds have now been 
placed on hadronic parityviolation in few-body reactions 
involving neutrons [Til . Il2j . With these developments, 
there is cause to study parity violation theoretically using 
QCD. Lattice gauge theory computations can determine 
the parity-violating couplings between hadrons from first 
principles. Indeed, the first exploration of hadronic 
parity violation using lattice QCD has appeared [liij . 
wherein a signal for the isovector parity-violating pion- 



nucleon coupling was obtained. While lattice QCD calcu- 
lations of parity violation are currently at quite an early 
stage, we are confident that precision information about 
hadronic parity violation will come from the lattice as 
refinements are made. To this end, one must understand 
the sources of parity violation at hadronic scales. The 
isovector channel has been studied at next-to-leading or- 
der [13]. In this note, we focus on the next-to-leading or- 
der corrections in the isotensor channel. In this channel, 
parity violation gives rise to couplings between two pi- 
ons and two nucleons, as well as parity odd pion-photon- 
nucleon interactions [15| . From a computational perspec- 
tive, this channel is expected to be statistically clean due 
to the lack of weak operator self-contractions. 

We present our findings as follows. First in Sec. [TT1 
we discuss isotensor parity violation at the scale of weak 
interactions, and determine the four-quark interaction 
in this channel at next-to-leading order accuracy in the 
QCD coupling. The renormalization group is then em- 
ployed in Sec. Mil to run the isotensor parity-violating 
interaction down to hadronic scales. As the evolution 
is renormalization scheme dependent, we present results 
for isotensor partiy violation using various renormaliza- 
tion schemes. We begin by obtaining results in 't Hooft- 
Vcltman dimensional regularization. Results for sev- 
eral regularization independent-momentum subtraction 
schemes are determined and compared in Sec. IIV1 which 
concludes this work. 



II. ISOTENSOR PARITY VIOLATION 

In the Standard Model, isotensor parity violation arises 
from the exchange of W and Z bosons between quarks. 
At scales below the W mass, My/, we can integrate out 
the weak vector bosons to arrive at an effective theory of 
the form 
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FIG. 1. One-loop QCD renormalization of weak-boson 
exchanges in the isotensor channel. The straight lines are 
quarks, while the curly lines are gluons, and the dashed lines 
represent W and Z bosons. Wave-function renormalization 
contributes but is not depicted. 



FIG. 2. One-loop QCD renormalization of the isotensor 
interaction in the parity- violating effective theory. Filled cir- 
cles denote the four-quark operator, while all other diagram 
elements are as in Fig. [T] Wave-function renormalization con- 
tributes but is not depicted. 



There is only one isotensor parity violating four-quark 
operator [l6j], namelj0 



O = (^r^^jLbrVgjt - ^(qr^q) L • {qT^q)i 



{l^r}. 
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The subscripts denote the left-handed and right-handed 
quark fields, qL,R — |(1 T 75)9, with the isodoublet 
q = (u 7 d) T . In the absence of electromagnetism, the op- 
erator O does not mix with other parity violating opera- 
tors. The corresponding Wilson coefficient, C(/x), evolves 
with the scale parameter fi so that the effective theory 
described by £p(r 2 is scale independent. 

At tree-level, the Wilson coefficient is scale indepen- 
dent and receives contributions from both W and Z 
exchange. To separate these contributions, we write 



q(o) _ qW _|_ Cyp, where the tree-level values of these 
coefficients are 



Cf = i-sin 2 . 
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At next-to-leading order, the Wilson coefficient of the 
operator O is determined by comparing the one-loop 
renormalization in the full and effective theories. For 
a complete discussion of QCD renormalization beyond 
leading logarithms, see [17] . The gluon radiation in 
the full theory is depicted in Fig. Q] After the quark 
wave-function renormalization has been accounted for, 



1 To minimize statistical noise in lattice calculations, the chiral ba- 
sis is not optimal because parity violation arises only upon taking 
the difference of hadronic matrix elements. A practical way to 
rewrite the isotensor parity- violating operator O for use in lattice 
calculations is O = — [(g r' A q) A (q r 3 q) v — ^(qfq) A ■ (qfq) v ], 
where (qq) a{qi)v denotes the Lorentz contraction of axial- vector 
and vector bilinears, (qq)A(qq)v = (<77n759)(97 A1 <?)- Written in 
this way, each term is parity violating. 



the diagrams shown are finite. In order that the effec- 
tive theory reproduces the full theory at one-loop, we 
must determine the matching coefficients by computing 
the effective theory diagrams shown in Fig. [2] These dia- 
grams require regularization, and we begin by utilizing 't 
Hooft-Veltman dimensional regularization with MS sub- 
traction. In this scheme, the axial-current vertex receives 
a finite renormalization at one loop [18], which indicates 
a non- vanishing two-loop anomalous dimension. To pre- 
serve chirality, we follow [l9l . [20j and augment the 't 
Hooft-Veltman scheme with an additional multiplicative 
renormalization chosen to force the two-loop anomalous 
dimension of the axial-vector current to vanish. Addi- 
tional renormalization schemes are considered in Sec. IIVI 
below. 

To perform the matching of the full and effective the- 
ories of isotensor parity violation, we choose to match 
at the scale [i = Mw- Diagrams with a W boson ex- 
changed can then be matched between theories without 
any explicit logarithms in the matching coefficients. To 
one-loop order, we find 



c w (m w ) = a 
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On the other hand, diagrams in the full theory with 
a Z boson exchanged will contribute logarithms of the 
form log (Mw/Mz) = log (cosf^) to the matching coef- 
ficients, and such contributions arise only from diagrams 
(c)-(/) in Fig. [T] Computing the matching coefficients 
for Z-exchange diagrams, we find 



Cz(M w ) = C 



(0) 
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Finally the desired initial value of the Wilson coefficient 
C(/z) at the VF-boson mass is simply 



C(M W ) = Cz(M w ) + C W (M W ), 
with the individual terms specified above. 



(6) 
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III. EVOLUTION 

The Wilson coefficient, C(/z), satisfies the renormaliza- 
tion group equation, 



^C(^)=r(M)C(A*), 



(7) 



where the anomalous dimension T(/i) has the expansion 
to second order 



a s f a A 2 
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with the scale dependence of T arising from the run- 
ning coupling, a s (/i). At next-to- leading order, the run- 
ning of the QCD coupling is determined by the equation 

- -a> m 2 - a (i#) 3 > ^th O = ii - §jv /5 

and /3i = 102- f TV/. 

The anomalous dimension T of the isotensor parity- 
violating operator, O, can be deduced from known re- 
sults concerning the QCD renormalization of |A5| = 1 
operators. Following [14j, we observe that parity invari- 
ance of QCD and of the renormalization scheme allows 
us to determine the anomalous dimension of O in Eq. ([2]) 
by considering the Vl <8> Vl combination of quark bilin- 
ears. Due to the isotensor nature of the operator, more- 
over, there are no penguin contractions to consider in the 
QCD renormalization. The anomalous dimensions aris- 
ing from current-current contractions of Vl ® Vl opera- 
tors have been determined in [21[ up to next-to-leading 
order. A closed basis of Vl ® Vl operators with only 
current-current contractions is 



o 2 = (V'iTm^M^'V^Ol, 



(9) 



where the ipi are each distinguishable quark fields, 
and the mixed brackets denote the color contraction 
(ipifa] [^3^4) = Sab ^1^2 V'lV'f- F° r the isotensor op- 
erator O, the analogous color rearranged operator O is 
identical to O after a Fierz transformation. Because 
the 't Hooft-Veltman scheme respects Fierz transforma- 
tions mixing of O into O is simply equivalent to 
a renormalization of O. With the anomalous dimension 
matrix Tq vll of the two operators in Eq. @, the desired 
anomalous dimension r of O is merely the sum of the 



first row of Fq vll . Hence we find 
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At next-to-leading order, the solution to the renormal- 
ization group equation is 



C{h')=U{ij!,h)C{h), 



with 
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where J = — Each of the quantities, /3o, Pi, and 
71, depends on the number of active flavors, Nf. This 
number changes as we cross heavy quark thresholds. As a 
result, the Wilson coefficient at hadronic scales, Aqcd < 
fi < m c , is given by 

C{n) = U {3) ([i,m c )U {i) {m c ,m b )U {5} (m b ,M w )C(M w ). 

There are no one-loop matching conditions required as we 
cross heavy quark thresholds. To this order, it is suffi- 
cient to account for these thresholds via the scale parame- 
ters A( N f) entering the definition of the two-loop running 

coupling, aa f \fj,,A^ N f^). The scale parameters are de- 
termined by enforcing continuity of the coupling at the 

heavy quark threshold, toq, namely a 



W (m Q) AW) = 



(mg,A< 



N, 



Consequently the Wilson coeffi- 



cient is described by a piece-wise continuous function, 
C(/x), of the renormalization scale fj,. 



IV. SUMMARY OF RESULTS 

Above we discuss parity violation in the isotensor chan- 
nel at the scale of weak interactions, and the QCD renor- 
malization group evolution to lower scales. Now we deter- 
mine the strength of the parity-violating isotensor oper- 
ator at hadronic scales in 't Hooft-Veltman dimensional 
regularization. Further schemes are discussed, and the 
corresponding results for isotensor parity violation are 
also determined in these schemes. 

With the value of the Wilson coefficient at the weak 
scale in 't Hooft-Veltman regularization, i.e. C{Mw) in 
Eq. ([5]), we can run down to hadronic scales using the 
solution to the renormalization group equation, Eq. (|13l) . 
This result is presented in Table HI and requires the MS 
masses of heavy quarks, and the value of a s {Mz), as 
described in pjj]. In this scheme, the next-to-leading 
order corrections are ~ 20% at the scale fx = 1 GeV. 

The conversion to other renormalization schemes can 
be achieved through one-loop matching. It is common 
to quote results in the naive dimensional regulariza- 
tion scheme, in which 75 remains anti-commuting in d- 
dimensions. At two- loop order, consistent results for the 
renormalization of four-quark operators can be obtained 
in this scheme, see, for example, [2l|. From one-loop 
matching, the Wilson coefficient of isotensor parity vi- 
olation using naive dimensional regularization is deter- 
mined in Table |U and differs by ~ 20% from that in 't 
Hooft-Veltman dimensional regularization. This is due to 
a relatively large matching coefficient compounded with 
the low scale \i = 1 GeV. The next-to-leading correc- 
tion in naive dimensional regularization is ~ 5% of the 
leading-order value. 

While dimensional regularization schemes are ideal for 
perturbative calculations, schemes that can be defined 
non-perturbatively are efficacious for lattice calcula- 
tions, such as the regularization independent-momentum 
scheme (RI/MOM) 22]. In this scheme, each external 



4 



TABLE I. Wilson coefficient of the |A7| = 2 parity- violating 
operator, C(/j,), at hadronic scales. We compare values ob- 
tained from leading-order (LO) evolution quoted in [15| to re- 
sults obtained in this work. At next-to-leading order (NLO), 
we compare results in various renormalization schemes. In 
each case, we scale the coefficient by its tree- level value, C^ ' . 
All momentum subtraction schemes employ Landau gauge. 



T O 


o v ) j c- 


LO [15] 


0.79 


LO 


0.70 


NLO 


C(lGeV)/C* (0) 


't Hooft-Veltman 


0.58 


Naive Dim. Reg. 


0.74 


RI/MOM 


0.77 


RI/SMOM( 7m ,^) 


0.67 


RI/SMOM( 7m , 7m ) 


0.75 


RI/SMOM((^) 


0.73 


RI/SMOM«, 7fl ) 


0.81 



leg of the four-quark operator is taken with an identi- 
cal off-shell momentum. Results in the RI/MOM scheme 
can similarly be determined from one-loop matching, and 
the resulting Wilson coefficient is given in the table. The 
choice of quark momenta in this scheme is quite excep- 
tional, for example, no momentum flows out of the opera- 
tor. It has been argued that an exceptional renormaliza- 
tion point is undesirable from chiral symmetry breaking 



considerations, and renormalization about an exceptional 
point generally exhibits poor convergence (23l . [2J|- To 
correct these maladies, schemes with a symmetric mo- 
mentum subtraction point (RI/SMOM) have been sug- 
gested in which momentum flows out of the operator. 
Results for the four different RI/SMOM schemes pro- 
posed for four-quark operators in (2f| (see also are 
given in Table HI 

Our study of isotensor parity violation shows that 
QCD evolution suppresses the magnitude of the inter- 
action at hadronic scales. Generally there is about 25% 
reduction, although this depends on the particular renor- 
malization. We compare several schemes, and find on av- 
erage ~ 10% differences arise, which gives an indication 
of the size of higher-order effects. It will be interesting 
to see if hadronic matrix elements are small in this chan- 
nel compared to the AI = 0, 1 channels. If this were the 
case, hadronic parity violation would parallel the AI = | 
rule, for which the maximal isospin changing channel is 
considerably suppressed. This is perhaps already sug- 
gested by the size of Wilson coefficients. Comparing re- 
sults in the 't Hooft-Veltman scheme from [lj], we see 
max{C A/=1 (lGeV)}/C A/=2 (lGeV) ~ 3. As there could 
be considerable cancellations between matrix elements in 
the AI = 1 case, we must ultimately wait for the evalu- 
ation of non-perturbative physics from lattice QCD. 
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